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Abstract 

This paper discusses the Space Shuttle Lift-Off CFD model developed for potential Lift-Off 
Debris transport for retum-to-flight. The Lift-Off portion of the flight is defined as the time 
starting with tanking of propellants until tower dear, approximately TO+6 seconds, where 
interactions with the launch pad cease. A CFD model containing the Space Shuttle and launch 
Pad geometry has been constructed and executed. Simplifications required in the construction of 
the model are presented and discussed. A body-fitted overset grid of up to 170 million grid points 
was developed which allowed positioning of the Vehicle relative to the Launch Pad over the first 
six seconds of Climb-Out. The CFD model works in conjunction with a debris particle transport 
model and a debris particle impact damage tolerance model. These models have been used to 
assess the interactions of the Space Shuttle plumes, the wind environment, and their interactions 
with each other and the Launch Pad and their ultimate effect on potential debris during Lift-Off. 


Introduction 

On January 16, 2003, the 28 th flight of the Shuttle Columbia was launched from Kennedy Space 
Center (KSC). This was the 113 th mission of the Shuttle Program, designated STS- 107. At 81.7 
sec after launch, when the Shuttle was at about 65,200 ft and traveling at Mach 2.46 (1,650 mph), 
a large piece of hand-crafted insulating foam came off die forward bipod area where the Orbiter 
attaches to the External Tank (ET). At 81.9 sec, the foam struck the leading edge of Columbia’s 
left wing 1 . There was a breach of the wing, eventually destroying the wing during reentry and 
causing the loss of Columbia and its crew on February 1, 2003. 

The need to determine the possible size, speed, impact location, and impact energy of debris led 
NASA and the Columbia Accident Investigation Board (CAIB) to use Computational Fluid 
Dynamics (CFD) techniques to analyze the event 2, 3 . The Space Shuttle Program also undertook 
the effort of identification and control of every possible source of debris liberation. The 
OVERFLOW CFD program 4 and chimera grid approach had been in use at NASA/JSC to 
simulate the flow over the integrated Space Shuttle Launch Vehicle (SSLV) for over 16 years. 
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These simulations were refined in support of the retum-to-flight effort. Through the leadership of 
the authors of Reference 2, significant geometric detail was added to the Vehicle CFD model 
such that the model consists of over 35 million grid cells. Steady-state CFD calculations were 
performed at several flight conditions duplicating the STS- 107 ascent trajectory. The CFD 
simulations provide not only the aerodynamic forces on the Vehicle but also a flow field for the 
ballistic integration of possible debris flight trajectories. 

In February 2004, as part of the retum-to-flight effort, NASA/MSFC was assigned the portion of 
SSLV flight regime from the beginning of propellant tanking up to the time when the Vehicle has 
cleared the launch tower. The active mechanisms for possible debris transport are gravity, 
Ground Winds, integrated propulsion system exhaust plumes and interaction of these mechanisms 
with Ground Facility structures. A new CFD model was constructed that included the pre- 
existing SSLV model from NASA/JSC augmented with all of the relevant Launch Pad Ground 
Facility structure. Point mass debris tracing was used for the aerodynamic transport of possible 
debris particles in this timeframe. This paper discusses the Space Shuttle Lift-Off CFD model 
developed for Lift-Off Debris Transport Analysis supporting retum-to-flight and continuing 
Space Shuttle operations. 


Debris Analysis and Risk Mitigation 


Aggressive debris mitigation and debris risk assessment have been carried out by NASA and its 
contractors supporting Space Shuttle retum-to-flight. Debris environment analysis for Lift-Off 
and Ascent follows a continuous improvement path to support all remaining Space Shuttle flights 
through International Space Station Assembly Complete. The Space Shuttle Program undertook 
mitigation of potential debris release and characterization of what debris might still be liberated 
on the basis of mass, material, density, shape, location and time of release. Potential debris items 
were compiled into a database and characterized as either expected or unexpected by design or 
process controls 5 . 

The most difficult requirement, derived from the greatest debris threat, for the Lift-Off Debris 
Transport Analysis to address is the ability to predict the velocity and direction debris might take 
if it is caught in an exhaust plume and turns back upward towards the SSLV. This is a 
probabilistic problem inasmuch as uncertainties in assumptions and models must be known and 
dispersions in the important transport variables must be understood. 

Foreign Object Debris (FOD) control at Shuttle Launch Complex 39 Pads A and B is a prime 
concern and a responsibility for all who work there. Corrosion from SRB exhaust products and 
the salt mist from the ocean incessantly weaken the structure. Rust flakes, bolts, nuts, and 
washers and other damaged hardware must be removed by continual Pad refurbishment and 
rigorous corrosion mitigation discipline. A disciplined mitigation effort was accomplished at Pad 
B before the retum-to-flight Missions STS-114 and 121 and is in progress at Pad A before 
launches resume there. FOD at the Pad is characterized as unexpected debris. 

The mitigation mass target limit for Facility rust and scale and, by inference, to any FOD on the 
Pad was 0.014 lbm, approximately the diameter and thickness of a U.S. 25 cent piece (a quarter). 
The MSFC Propulsion Systems Engineering and Integration Office is responsible for 
determination of the lift-off debris environment. MSFC Engineering developed the lift-off CFD 
model and is performing the lift-off DTA. 
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The Lift-Off CFD model requirements are extensive. The model shall be applicable to the 
Tanking period prior to launch, while ice or other potential debris might fall on the Vehicle. It 
must be applicable to Main Engine Start, at Solid Rocket Booster Ignition, and during Climb-out 
past ‘Tower Clear’, where vibration and plume blast may result in the liberation of debris that can 
impact the Vehicle. 


The Nominal Launch Sequence 


Shuttle launch operations are to be simulated beginning with the tanking period when cryogenic 
propellants are on-board and dropped to the SSMEs for chill-down and recirculation. Ice 
formation, which depends on temperature, humidity and ground winds, presents the threat of 
falling ice debris from this time forward. Ice may come from the External Tank (ET), feedline 
bellows and brackets, aft fittings, and main engines. 

At ~ TO - 6 sec, the Space Shuttle Main Engines (SSMEs) are started. By ~ TO — 2 sec the 
SSMEs have reached 100 % rated power level. SRB Ignition Command is issued at TO at which 
time the SRB holddown explosive fasteners are fired and the TO umbilicals are retracted. The TO 
umbilicals consist of the L02 and LH2 Tail Service Mast (TSM) umbilicals on the Mobile 
Launch Platform (MLP) and the GH2 Vent Line at the ET Intertank. Tne time ~ TO + 0.3 sec is 
the moment that the SSLV begins to rise as it is the time that the SRBs have attained full thrust. 
During the climb-out period, the SSLV is under the threat of ‘Pad debris’ that may be blown 
upward as a result of plume interactions with Ground structures. The SSLV is under the threat of 
falling debris that can fall from either the SSLV or Pad structures above. The SSLV vibration, 
acoustics, and ignition overpressure may excite debris on either the SSLV or the Pad to be 
liberated at any time after SSME Start. 

There is therefore a time during SSLV Climb-out after which the interactions with Launch 
Facility structures no longer can have an influence on the SSLV debris environment. Upward- 
traveling debris particles can no longer overtake the SSLV, and falling particles from the highest 
locations on the Launch Access Tower can no longer have a possible trajectory to the SSLV. We 
define this ending time for lift-off DTA to be ‘Tower Clear’. It is approximately when the SRB 
aft skirt aft extremity is higher than the top of the Fixed Service Structure (FSS) Access Tower 
and is before ~ TO + 6 sec. 

We chose to perform CFD simulations using steady state time slices, referred to as Cardinal 
Points (CardPts). The Cardinal Points are shown in Table I and were chosen based on knowledge 
of key events in the time sequence. The first, CardPt 1, can be any time when propellants are 
loaded and the Vehicle and SSMEs are fully chilled and there is ice/frost formation. CardPt 2 is 
after the three SSMEs have attained 100 % rated thrust level. CardPt 3 is reserved for future 
analyses of the SRB ignition transient (not a part of the present steady-state simulations). 
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Table I. CFD Cardinal Points 


Wind Speed 

Time 

6.7 knots 

20.1 knots 

33.5 knots 

Wind Direction 

sec from TO 

45 

90 

135 

180 

225 

270 

315 

360 

45 

90 

135 

180 

225 

270 

315 

360 

45 

90 

135 

180 

225 

270 

315 

360 

CardPtl 

prior to -6.0 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

CardPt2 

-2.0 

X 

X 




X 

X 

X 

X 

X 




X 

X 

X 

X 

X 




X 

X 

X 

CardPt3 

- 

Reserved for Later Use 

CardPt4 

0.3 

X 

X 




X 

X 

X 

X 

X 




X 

X 

X 

X 

X 




X 

X 

X 

CardPt5 

1.9 








X 

















CardPt6 

3.0 








X 

















CardPt6a 

3.0 








1 

















CardPt7 

4.0 








X 

















CardPt8 

5.0 








X 

















CardPt9 

6.0 








X 

















CardPt9a 

6.0 
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Note 1 : Cardinal Points 6a and 9a used forward flight velocity from the STS-1 1 3 BET. 


CardPt 4 is at the moment of lift-off, at ~ TO + 0.3 sec. CardPt 5 is at a critical time in the plume 
flow transient when the three SSME plumes are fully contained in the SSME exhaust hole in the 
MLP but the two SRB exhaust plumes begin to impinge on MLP structure, at ~ TO + 1.9 sec. As 
the Vehicle climbs, the SRB plumes begin to impinge on the SRB holddown posts and haunches 
starting at ~ TO + 1.5 sec. This SRB plume impingement continues on but begins to abate by ~ 
TO + 4 sec as the Vehicle gains more height. The impingement pattern in the SRB exhaust holes 
changes during the 1.5-4 sec interval because the Vehicle is rising and also because it drifts in 
the Vehicle body axis plunging direction, Pad north (Vehicle -Z T coordinate direction), due to the 
thrust vector of the SSMEs. 

A new flow phenomenon is initiated at ~ TO + 3 sec as the SRB plumes are larger than the width 
of the SRB exhaust holes and begin to spill out onto the MLP Deck. Gross SRB plume spillage 
and impingement on the flat Deck continues on until well after ‘Tower Clear’. This flow 
condition is captured in CardPt 6. 

At ~ TO + 4 sec, there is considerable interaction among the five exhaust plumes and the Facility 
structures, but the plumes from the two lower SSMEs begin to impinge on the TSMs and spill out 
of the SSME exhaust hole. This flow condition is captured in CardPt 7. 

Simply continuing with the climb-out phase on past ‘Tower Clear’ at even 1 sec intervals, 
CardPts 8 and 9 have the Plume-Facility interactions at TO + 5 sec and 6 sec, respectively, and 
‘Tower Clear’ occurs at ~ TO + 5.5 sec. 


Lift-Off Debris Transport Analysis 


The Lift-Off Debris Transport Analysis process is illustrated in Figure 1. The analysis begins 
with a CFD simulation of the flowfield about the Vehicle and related Ground Structure. The 
prediction is at a given time pre-launch, during lift-off or ascent flight and is a quasi steady-flow 
solution representing an instant in time. The time slice for the CFD analysis is referred to as a 
Cardinal Point (CardPt). 

Initial debris locations and release velocity vectors are known or estimated. The debris is 
modeled as a point mass where just its material and density, mass m, drag coefficient, Co, cross- 
section reference area, A ref , and therefore its ballistic number, BN are required. A debris transport 
trajectory integration program, debris 6 , computes the trajectory of the debris particle through the 
CFD flowfield from the point of release either to impact or miss the SSLV. 
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Impacts are detected and characterized as to their impact relative velocity, kinetic energy, and 
angle by the dprox 1 program. Many debris particles of different size and shape may be released 
as point masses into a given CFD flowfield and the maxima in impact kinetic energy or angle or 
velocity sorted in a given impact zone on the SSLV. Many such sorts for the impact quantity 
maxima make up the potential SSLV debris impact environment. 



Figure 1. The Debris Transport Analysis Process 


AnalysisTools 

The OVERFLOW CFD program is used to simulate the Lift-Off flowfield because the debris and 
dprox tools were developed to use the CFD flowfield from this particular CFD solver. Having a 
full existing tool suite and a common approach for Lift-Off and Ascent Debris Transport Analysis 
made it a pragmatic decision to use the OVERFLOW program. The time required to complete 
the reconstruction of an equivalent tool suite using based on another CFD program was judged to 
be prohibitive. The only risk associated with this decision is that OVERFLOW uses a ‘density- 
based’ algorithm and had only recently been upgraded to include a preconditioned formulation 
for low-Mach number accuracy. It was unknown whether this formulation was capable of 
robustly and accurately simulating this problem which depends on robust and accurate flow-field 
predictions with significant portions of the same simulation having Mach numbers ranging from 
0.01 to 6. 

Previous CFD models 8, 9 of the SSLV On-the-Pad have been only Pre-launch. They were of 2D 
section cuts through the Vehicle and of an External Tank (ET) - Solid Rocket Booster (SRB) 3- D 
half-model with pre-launch winds and humidity influences on the ET chill effect and utilized the 
Rockwell International Unified Solution Algorithms code 10 . There have been wind tunnel tests 
measuring Vehicle integrated force and moment coefficients from pressure-area loads due to pre- 


5 



launch Ground Winds 11 . This and prior 12, 13 wind tunnel data are being used to support validation 
of the lift-off CFD model. 


Objective and Approach 


The objective was to develop a Lift-Off CFD model in which the transport of potential lift-off 
debris could be evaluated. This was in order to provide impact data from which debris damage 
tolerance assessments could be made as to location and severity of debris impacts. Rapid 
response was needed to support return-to-flight. There were several important simplifications 
made in order to be responsive. 


Simplifications 

The simplifications are rank-ordered in what is believed to be the order of their importance to 
Lift-Off DTA, and priorities for continuous improvements were set based on this rank order. 

The first simplification was to use the quasi steady state approximation. 

It is not yet possible to simulate using CFD the SSLV lift-off transient in a time-accurate mode 
and include the necessary geometric detail. The SSLV moves with respect to the Launch Pad, 
various components articulate and power setting changes occur with the SSLV in flight as a 
function of time. Therefore we resorted to quasi-steady simulation using the CardPts and Ground 
Wind dispersions within a given CardPt. Dispersions cover an envelope of possible Ground 
Winds day to day (see Table 1). This decision largely defined the required number of CFD 
simulations which wound up being sixty-one. 

The time scale between the release of a piece of debris and its impact on the SSLV is on the order 
of milliseconds when the Vehicle is flying at supersonic speeds (as in the case of the STS-107 
foam impact). The quasi steady-state approximation for the CFD flow field is justified in that 
case. However, the time for debris being released at lift-off, being accelerated by gravity and/or 
interaction with the propulsive flows and ground facility, and coming back to an impact on the 
SSLV may be on the order of seconds. The SSLV will have moved with respect to the Facility 
during that elapsed time from the debris release to its impact. For Lift-Off debris transport 
analysis, the quasi steady-state approximation for the flow field is clearly not justified. 

We propose that by selecting a large enough number of time slices (CardPts) it is possible to 
represent the real transient problem. This is an approach that may bound the debris transport 
phenomena. Comparing CFD results with launch imagery from many flights has been the 
verification approach taken to try to verify that we have enough CardPts and that we have 
included sufficient geometric detail. But this approach is no substitute for striving to work the 
transient problem in the future based on our continuing improvement philosophy stated above. 

The second simplification was to use a single fluid specie, air as a perfect gas, with y= 1 .4. This 
is the working fluid everywhere; hot air at high velocities is used to simulate all exhaust plumes. 
Thus, the CFD model approximates the plume gaseous composition and neglects completely the 
condensed-species particulates in the SRB exhaust plumes. State-of-the-art SRB plume 
simulations 14 (single plume, axisymmetric) include the particulate flow, but two-phase flow 
simulation is a complexity that would taken too long to implement in this 3-D model. This 
simplification still provided reasonable plume jet velocities and pitot impingement pressures, but 
left out ~ 40 % of the SRB exhaust mass flow. 
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Real gas plume simulations have been evaluated for Space Shuttle ascent flight conditions 15 using 
OVERFLOW with the conclusion that y= 1.4 plumes with MW = 29.0 would provide adequate 
aerodynamic simulation with the plume expansion (and blockage effect) at altitude based on the 
Orbiter surface c p . 

For the lift-off regime, accurate simulation of shocks, expansions, and mixing have considerably 
different modeling requirements 16 than do plumes at altitude which are underexpanded. Real gas 
thermal and fluid transport properties must be preserved to retain all of the plume energy properly 
and capture the aspiration, impingements and reverse flows accurately. In the spirit of continuing 
improvement, this simplification should also be amended. 

The third simplification was to de-feature some of the Facility structure surface details. For 
example, the Fixed and Rotating Service Structures were neglected. The SRB Exhaust Holes in 
the Mobile Launch Platform (MLP) were significantly de-featured, but retain the basic 
geometrical shape of the actual hardware. What results is considered to be a first-order geometry 
representation for effects of SRB plume impingement and recirculation was constructed. Even 
with aggressive de-featuring to manage the size of the CFD grid, some CardPts contained up to 
170 million cells. 

Unstructured grid technology offers the flexibility to cluster more grid points where needed and 
less than where they are not needed. Converting to a CFD solver which supported unstructured 
grids would allow the addition of more ground structure features in a continuous improvement 
path. 

The fourth simplification was to neglect the Sound and Ignition Overpressure suppression water 
bags and sprays. This decision eliminated > 1,000,000 gpm sound suppression water and steam 
flow from the total flow. The water sprays are below the MLP elevation considered most 
important for debris sources at lift-off and post-liftoff, and a blocking filter was used in debris to 
restrict any debris sources coming from below the water sprays locations. Post-liftoff sprays 
come on for noise suppression at the top Deck of the MLP after ~T0 + 3 sec. The possible effect 
that energy extraction from the plumes and the volume expansion from flashing of the water to 
steam could have on plume entrainment and/or rebound of debris particles was not included 
because it was judged to be small, 

The fifth simplification was a simple slug flow profile Ground Wind that did not have the Earth 
boundary layer effects of vertical wind shear velocity profile and inclusion of the berm effect of 
the ‘Hill’ close to the Pad. Ground winds were simply applied as slug flow dispersions every 45 
deg around the compass, at three wind speeds - mean, near 2<xhigh, and near 3 <7 high. 

This simplification left out the Ground Wind blockage effect of the Fixed Service Structure. This 
blockage will have an influence on debris transport particle trajectories when the Ground Wind 
direction is near either 270 deg (from the West) or 90 deg (from the East), but was considered to 
be small based on the wind tunnel results of Reference 1 1. 


Grid Construction. Vehicle Configuration, and Boundary Conditions 

The present model combines the Space Shuttle body-fitted structured grid with a newly 
constructed body-fitted grid for the MLP and Ground structures. The SSLV model has all 
protuberances on the ET resolved down to y + = 1 and over 80 million grid points including 
exhaust plume grids. It corresponds to the STS-1 14 configuration and has 270 multi-block grid 
blocks. Plume boundary (‘ start line’) conditions are applied at the nozzle exits of the SSMEs and 
SRBs and are obtained from the JANNAF Standard Plume codes 14 . The SSME nozzle flow 
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definition agrees well with solutions obtained using other codes such as the Finite Difference 
Navier-Stokes (FDNS) solver 16 , for instance. 

A detailed Launch Complex 39 Facility CAD model provided the geometry needed for generating 
the lift-off CFD grid model 17 . Integrated with the Launch Facility, the CFD volume grid was as 
large as 170 million grid points. 


Low Mach Number Pre-conditioning 

Pre-conditioning was necessary for OVERFLOW code operation at Ground Wind velocities. 
There is the extreme range in velocities all the way from nozzle exit plume velocities (Mach 6) 
down to the Ground Wind free stream ambient velocities (Mach 0.01). A separate study 18 was 
conducted at the start of the Lift-Off DTA effort to verify, validate, and test the low Mach 
number pre-conditioning extension implemented in the OVERFLOW 2.0s flow solver. The study 
concluded it was possible to demonstrate good performance of the code down to 0.01 Mach 
number (ambient wind speed of 6.7 knots on a standard temperature-pressure day) for low-speed 
steady-flow simulations. This was with the Spalart and Allmaras turbulence model 19 at Re up to 
2.5 X 10 7 . 


The Cardinal Point Cases 

Sixty-one individual CFD cases were executed at the operating conditions for nine CardPts, Table 
I. CardPts were selected pre-launch just before SSME Start and then basically at 1 sec intervals 
at ~T0 - 1 sec, -TO + 0.3 sec, and TO + 2, 3, 4, 5, and 6 sec. There were Ground Wind 
dispersions at 24 wind speed-direction combinations pre-launch. Three wind speeds were 
selected - 6.7 kts, 20.1 kts, and 33.5 kts (corresponding to M in f= 0.01, 0.03, and 0.05). Wind 
directions selected were every 45 deg around the compass. See Table I for these details. 

SSLV flight speed, altitude, and drift, rocket chamber pressures (Pc’ s), gimbals and elevon 
deflections were taken from a typical mission to the International Space Station flight trajectory 
reconstruction. Two CardPts left out the Ground structures in order to include the Vehicle 
forward flight velocity, CardPt 6a at TO + 3 sec and CardPt 9a at TO + 6 sec. The CFD model 
included articulating SSME gimbal angles and elevon surface deflections as they are commanded 
in flight. These changes in SSLV geometry and times at which they occur are listed for each 
CardPt in Table II. These schedules were followed closely, not precisely, in the CFD model. 


Table II. SSME Gimbal and Grbiter Elevon Deflection Schedules 



Time 

sec from TO 


SSME 2 and 3 
Gimbal, deg 

Inboard Elevon 
Defl., deg 

Outboard Elevon 
Defl., deg 

ram* 

- - 8 

0 yaw, 16 pitch 

0 yaw, 10 pitch 

0 

0 


-2.0 

0 yaw, 16 pitch 

0 yaw, 10 pitch 

0 

0 

BESM 

- 

0 yaw, 16 pitch 

0 yaw, 10 pitch 

0 

0 

mssm 

0.3 

0 yaw, 16 pitch 

0 yaw, 10 pitch 

0 

0 

BB 

1.9 

0 yaw, 16 pitch 

0 yaw, 10 pitch 

0 

0 

ran* 

3.0 

0 yaw, 15 pitch 

0 yaw, 9 pitch 

1 

1 

fSSffl 

3.0 

0 yaw, 15 pitch 

0 yaw, 9 pitch 

1 

1 

BMS1II 

4.0 

0 yaw, 12.5 pitch 

0 yaw, 7.5 pitch 

3 

3 

KMHI 

5.0 

0 yaw, 12.5 pitch 

0 yaw, 7.5 pitch 

5 

5 


6.0 

0 yaw, 12.5 pitch 

0 yaw, 7.5 pitch 

7 

7 

IBSBi 

6.0 

0 yaw, 12.5 pitch 

0 yaw, 7.5 pitch 

7 

7 
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Between TO + 2.5 sec and TO + 7.5 sec, the Orbiter elevens are commanded to slew at 2 deg/sec 
from a preset 0 deg to 9 deg down on the outboard elevens and 10 deg down on the inboard 
elevons. At TO + 2 sec the SSMEs, initially set yaw parallel and at 16 deg up for the upper 
engine (Engine 1) and 10 deg up for the two lower engines (Engines 2 and 3), are commanded to 
pitch down. By TO + 4 sec, the upper engine is pitched down 3.5 deg to a position of 12.5 deg, 
and the two lower engines are pitched down 2.5 deg to the 7.5 deg position. 


Grid Sequencing and Convergence Criteria 

Grid sequencing was used to reduce the computational effort required to arrive at a steady state 
solution. This was accomplished in a seven-step run recipe wherein the flow around the SSLV 
and Ground is developed over 1000 iterations without the plumes activated. The second step uses 
the ‘slow-start’ feature within OVERLFOW to ramp up the plume boundary condition over 500 
iterations and then develop the plumes with an additional 1000 iterations. The first two steps are 
accomplished on a 4 th level sequence grid. Then the solution is refined in step 3 three by running 
1000 iterations on sequencing level 3 and so on until the solution is proceeding at 2 nd order at grid 
level 1 on all grids except the SSME grid box, which could be run only at 1 st order. Depending 
on CardPt and wind speed, the total number of iterations required for convergence ranged from 
20,000 to 1 10,000. Grid sequence level 1 is the actual fine-mesh grid. Grid sequence level 2 is 2 
X coarser, grid level 3 is 4 X coarser, and grid level 4 is 8 X coarser than the fine grid. 

In addition to watching the residuals reduce five orders of magnitude or more, the force and 
moment coefficient (FOMOCO) integration utilities in OVERFLOW were exercised to watch the 
integrated Shuttle Vehicle, each Element, and selected structural component FOMOCOs settle to 
constant values with noise having decayed to a minimum. 

The Space Shuttle Vehicle and Launch Facility CFD solutions were performed using a 
combination of NASA/MSFC Linux clusters and the Columbia SGI Altix system at NASA/ARC. 
The debris tracing and impact detection and characterization programs were executed on a 
NASA/MSFC SGI Origin 3900. A full run according to the recipe above could be accomplished 
in 6 days using 128 epu’s of the Columbia supercomputer. Ground wind dispersion cases were 
usually started from the solution to another ground wind direction and required about 21 days of 
32 epu’s on an Athlon-MP based Linux cluster. 


CFD Flowfield Results 


For the first time, a full 3-D CFD model has been developed for the Shuttle before and during 
Lift-Off. The SSLV and Ground Facility models have been combined and the Ground Winds 
have been modeled. The SSLV is moved in a series of steady state simulations above the Ground 
to simulate the plume flow interactions with Facility structures and what effects those interactions 
may have on potential lift-off debris transport to impact on the SSLV. 

In Figure 2 which corresponds to CardPt 1, the wind originates from the East and wraps around 
the vehicle, causing local areas of flow separation and recirculation which could affect debris 
trajectories. The downstream side of the SSLV is almost completely separated flow as shown by 
the velocity vectors and synthetic oil flow patterns in Figure 2. Two independent sets of 
experimental data 11-13 exist for the wind-induced forces on the SSLV on the launch Pad. In 
Figure 3, the data sets best corresponding to the CFD model geometry, i.e., without the Fixed and 
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Rotating Service Structures (ESS and RSS, respectively) are compared to the normal and side 
force coefficients computed from CardPt 1 from all wind speeds and directions. The computed 
force coefficients reasonably replicate the experimental data in both value and nature of response 
to different wind directions. While this result falls far short of validating the predicted wind 
interactions with the SSLV and ground geometry, it does represent one important component of a 
rigorous validation argument. 

Figure 4 shows the predicted wind patterns for CartPt 2, the case where only the SSMEs are 
operating. In contrast to CardPt 1 flow patterns, the presence of the aspiration from the SSME 
plumes radically reduces the separated flow volume. Velocity vectors in Figure 4 indicate a 
general downward flow direction imparted to the flow localized around the SSLV. Synthetic oil 
flow patterns illustrate the influence of the SSME plume aspiration near the SSLV body. The 
flow patterns converge on the SSME plumes and could influence debris trajectories towards the 
base of the SSLV as compared to their trajectories in a CardPt 1 flow field. 

To illustrate the consequences of the Second Simplification above, Figure 5 compares a single 
SSME plume in a near quiescent flow-field computed using finite rate chemistry from JANNAF 
standard codes to that computed using the OVERFLOW CFD program with the second 
simplification active. The OVERFLOW computation was performed using settings that rendered 
the algorithm to be first order spatially accurate. This was necessary because the OVERFLOW 
solver produced IEEE math exceptions when the second order algorithm was used. Significant 
interactions with the developers of the OVERFLOW program over a significant period of time 
did not result in a solution to this issue. The net result is that the prediction of the SSME plume 
used in this model degrades with increasing distance from the SSME nozzle exit plane. Figure 5 
suggests that the SSME plume prediction is reasonably accurate only to approximately ten feet 
downstream of the nozzle exit plane. Clearly, this is an aspect of the CFD model that needs to be 
improved. 

Figure 6 shows the predicted wind patterns for CartPt 4, the case where all 5 plumes are 
activated. In contrast to CardPt 2, the addition of the SRB plumes further strengthens the 
downward flow and almost completely collapses the separated regions of the flow around the 
SSLV. The synthetic oil flow patterns indicate a strong downward flow over most of the near 
body surface, converging towards the plume boundary layer surfaces. These flows towards the 
vehicle can influence debris trajectories to impact the SSLV both more often and with more 
energy than gravity and ground winds alone. 

A comparison of the SRB plume computed from standard JANNAF codes is compared to that 
computed using OVERFLOW with the second simplification active in Figure 7. For the first 50 
feet downstream of the nozzle exit, the two flow fields agree reasonably, but not perfectly. In 
contrast to the SSME plumes, the OVERFLOW program is capable of robustly simulating the 
SRB plume using a second order accurate spatial algorithm. Thus the differences between the 
two plumes in Figure 7 are the appropriate fluid species properties and the lack of the particulate 
model in the OVERFLOW simulation. Both simplifications must be addressed in the spirit of 
continuous improvement. 

Figure 8 illustrates the impact of the plume aspiration on the flow-field, in part a) during CardPt 2 
with the just the SSMEs active and in part b), with all five plumes active. Isosurfaces of constant 
Mach number around the base of the SSLV illustrate the acceleration and convergence of the 
flow. Any debris liberated in this region can be exposed to plume aspiration velocities on the 
order ofM = 0.1 or 115 feet per second accelerating the debris towards the Vehicle itself. 
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When the Vehicle has risen 20-30 ft or so above the Pad, the Vehicle’s drift to the north causes an 
asymmetric SRB plume impingement and recirculation pattern to set in the SRB Exhaust Hole 
below. CardPt 5 is definitely into this flow condition with the SSLV height at ~ TO + 1.9 sec 
being ~ 27 ft and the drift being ~ 3-4 ft. Figure 9 shows some close-up views of the SRB 
exhaust plume flow from CardPt 5. On the upper left is the streamline flow close in around the 
SRB plumes. The streamlines go into the SRB Exhaust Holes (aspiration). The plume isotherm 
contour, colored by Mach number up to Mach 1, shows approximately the boundary of the plume 
supersonic core. The view on the upper right shows details of the SRB holddown structure inside 
the SRB Exhaust Hole. The surfaces of the four holddown posts and support haunches are 
colored by the computed plume impingement pressure. The impingement pressure is 23 psi or 
more on the north holddown posts and haunches (red) as a result of the SSLV northerly drift. 

The flat surfaces of the haunches below the SSLV provide particle recirculation and rebound 
surfaces and the side walls provide possible surfaces for particle ricochet back upwards toward 
the SSLV. This flow condition, which begins at ~ TO + 1.5 sec and persists well past ~ TO + 3 
sec, is therefore a critical time for the potential upward transport of liberated debris particles. The 
SRB plume impingement effect causes upward flow pockets with as great as 650 fps upward 
velocity maximum found (colored regions > 50 fps illustrated) around the south comers of the 
SRB Exhaust Holes and on the TSM north faces, beneath the Orbiter, and at the center of MLP 
Deck 0, beneath the ET. These pockets are illustrated in the bottom right comer of Figure 9. The 
image on the bottom left hand comer of Figure 9 shows that the SRB plumes have not yet begun 
to spill out of the SRB Holes on the MLP deck. 

There is a change in the nature of the flow field beginning at ~ TO + 3 sec. The SRB plumes 
grow in diameter with increasing distance from the nozzle exit plane and they begin to spill out of 
the SRB Exhaust Holes and onto the MLP Deck. This flow condition is illustrated in Figure 10, 
which is the flowfield at ~ TO + 3 sec. The SSLV height is ~ 67 ft and the drift ~ 7 ft at this time 
and the SSLV has attained ~ 48 fps forward flight velocity and corresponds to CardPt 6. 

At the upper left in Figure 10 is the same plume flow plus streamlines view as before, but the 
streamlines spill outward around the plume periphery onto the MLP deck. There are significant 
horizontal components of the flow out across the deck as seen in the view from above at the lower 
left. On the upper right, the plume impingement pressure (red) covers a larger surface area on the 
holddown posts and haunches and will lessen as the SSLV continues to climb. 

The pockets of upward flow > 50 fjps are smaller and it can be seen that the plume spillage causes 
horizontal vectors causing the flow recirculation to be directed more outward and away from the 
Vehicle direction. This is a condition of expected lessening severity for debris particle upward 
transport as the velocity field now tends to turn particles liberated after ~ TO + 3 sec outward and 
away and the Vehicle is higher and accelerating away from the threat of upward-traveling debris. 

CardPt 6a flowfield views are shown in Figure 11. Plume flow cross-section cuts are shown in 
the upper left and upper right views. The color contours are Mach number in the plumes, highest 
in the SSME plumes. The upper view in the center shows mass-less streamlines colored by 
density going from free-stream (red) to the SSME plume boundaries (blues) and indicating the 
very strong aspiration into the orbiter base region. The indication is that debris particles drawn 
toward the base will tend to be drawn to highest velocities at the body flap, elevons, and the 
SSME nozzle bells. 

The view at the lower left in Figure 1 1 shows a typical set of liberated debris particle hit locations 
from the debris and dprox codes run with the CardPt 6 CFD flowfield. The debris particles may 
have come from above or from below to result in impacts at these locations on the Orbiter. 
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The view at the lower right in Figure 11 shows velocity contours in a transverse cut through the 
SSME plumes and just above the plane of the SRB nozzle exit planes. If upward-traveling debris 
makes it past the blockage volume of the SRB exhaust plumes, it still has to negotiate the strong 
aspiration down-flow of the SSME plume aspiration to reach the Orbiter aft extremity regions 
from below. 

These CFD results illustrate flow characteristics known to exist at launch based on the fact that 
they show features of the flow seen in launch film imagery and evident in bum and damage 
patterns left on the Facility structure after launch. 


Cycle-2 CFD Flowfield Model 


The strong interactions of the SRB plumes with the SRB holddown posts, resulting in upward 
flow primarily motivated a follow on effort to refine the fidelity of the CFD model. This effort 
has been termed the Cycle-2 Lift-off Debris Transport Analysis effort and has been underway 
since the beginning of fiscal year 2006. The Cycle-2 effort includes improvements in all aspects 
of the debris transport analysis; the CFD-model related improvements will be briefly discussed 
here. 

It is still necessary to make simplification #1 in the Cycle-2 effort. However, Simplification 2, 
the use of air as the plume effluent, can be improved drastically. We have chosen to use the Loci- 
CHEM CFD program 20, 21 , one of the production-mode CFD tools at NASA/MSFC. This choice 
allows the SSME plume to simulated with both a quite realistic effluent and using a second order 
spatially accurate algorithm. Eulerian and Lagrangian-based particle models are being 
implemented primarily to model the SRB internal and plume flows and will be available soon for 
use in the Cycle-2 effort. 

It is Simplification 3, the geometry simplification, which is receiving the most attention. The 
Loci-CHEM CFD program uses unstructured grids, which makes the construction of grids for 
complicated, detailed geometry much easier. Figure 12 shows the detail of the SRB Hole and 
Holddown posts available in the CAD mode, the Cycle- 1 CFD model and the Cycle-2 CFD 
model. The SRB-Holddown post interactions will be investigated by decoupling the SRB plume 
and SRB ground model from the rest of the SSLV and ground facility. In this manner the upward 
flow predicted from the Cycle- 1 model in Figures 9 and 10 will be assessed using significantly 
higher fidelity physical and geometrical models. 

Recently, the Loci-CHEM program has been upgraded to support real fluid equations of state, 
which allows the modeling of the Ignition Overpressure and Sound Suppression Water flows (see 
e.g. References 22 and 23, The use of this new capability will be attempted as part of the Cycle-2 
effort to address Simplification 4. 


Summary and Conclusions 


For the first time, a comprehensive CFD model of the SSLV during Lift-Off has been constructed 
and executed. This CFD model contains a detailed model of the SSLV and couples it with a 
detailed model of the MLP, Flame Deflector, and Exhaust Trench, and Ground geometry. In 
addition, the effect of the wind environment has been captured by evaluating three wind 
velocities every 45 degrees around the compass. The wind environment and plume aspiration 
effects have been coupled to arrive at a prediction of the flow field through which potential debris 
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trajectories are integrated. Over 60 CFD simulations have been completed in Cycle 1 with 
hundreds of debris simulation cases run with tens of thousands of potential debris transport 
trajectories generated. The knowledge of the nature of the interactions of the plume-generated 
flow field with the Ground geometry has been used to further focus potential FOD removal and 
mitigation. In particular, the identification of high speed upward flow, capable of driving debris 
upwards towards the SSLV with significant velocities, has been identified. The Cycle 1 debris 
database has been built for expected debris and is being used for Day of Launch real-time mission 
support. A basis for real time mission decision support has been gained as a result of knowledge 
from the lift-off CFD flowfield analysis. 

The Cycle-2 effort, largely motivated by the potential for driving debris upwards towards the 
SSLV, is now underway. Cycle-2 systematically addresses the Simplifications required for the 
construction of the Cycle- 1 CFD model. In particular, the fidelity of the plume simulation details 
and the fidelity of the Ground geometry from which the interaction of the SRB plumes cause the 
upward flow have been increased. 
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Nomenclature 


BET Best Estimated Trajectory (of the Vehicle at Lift-off, reconstructed from flight data) 
BN Ballistic number (of a debris particle), = m/Cj>4 re / 

CFD Computational Fluid Dynamics 
Cd Drag coefficient (of a debris particle) 
c p Pressure coefficient = (pi 0C J -pmj)/q,nf 
DTA Debris Transport Analysis 
ET External Tank 
FOD Foreign object debris 
FSS Fixed Service Structure 

Y Ratio of specific of heats of gas species (ideal gas assumption) 

JANNAF Joint Army Navy Air Force 

IEEE Institute of Electrical and Electronics Engineering 

MLP Mobile Launch Platform 

MW Molecular weight (of plume gas species) 

PAL Protuberance Air Load (on die ET) 

Pc Chamber pressure (of the SSME and the RSRM), psia 

p Static pressure, psia 

q Dynamic pressure, psi 

RSRM Reusable Solid Rocket Motor 

Re Reynolds number 

SRB Solid Rocket Booster 

SSLV Space Shuttle Launch Vehicle 

SSME Space Shuttle Main Engine 

<j Standard deviation 

TO Time of SRB Ignition Command 

TSM Tail Service Mast 

t Time (in a given launch trajectory as measured from TO) 

X T Shuttle Tank axial coordinate (positive going aft) 

Y t Shuttle Tank lateral coordinate (positive out RH Wing) 

Z T Shuttle Tank vertical coordinate (positive out Vertical Tail) 
y + Non-dimensional height in the turbulent wall boundary layer on a surface 
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Subscripts 

i Sequence of Cardinal Point Cases (given Vehicle trajectory, given Ground Wind 
condition) 

j Particular flight traj ectory (V ehicle climb-out) 

inf Conditions at infinity (far field) 
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Figure 2. Flow field predicted from CardPt 1 with wind speed of 6.7 knots and wind direction of 90 

degrees (from the East). 
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Figure 4. Flow field predicted from CardPt 2 with wind speed of 6.7 knots and wind direction of 90 

degrees (from the East). 
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Figure 5.Comparison of single SSME plume computed using finite rate chemistry (left) and computed 

using OVERFLOW and air as the working fluid. 
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Figure 6. Flow field predicted from CardPt 4 with wind speed of 6.7 knots and wind direction of 90 

degrees (from the East). 
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Figure 7 Comparison of SRB plume computed using standard JANNAF codes and the OVERFLOW CFD program. 
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Figure 8. The Potential Flow Nature of the Flow at Lift-Off 
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Figure 9. Cardinal Point 5. 
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Figure 10. Cardinal Point 6. 
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Figure 11. Vehicle Climb-out Solutions in the Timeframe after ~ TO + 3 sec 






Figure 12. Improvements to SRB Exhaust Hole Geometry 




NASA 



"D 

CO 

Q. 


CO 


Q) (/) 

■o w 

I - 
« (0 

o < 

E * 

(0 o 

c 9- 

(0 


s? 


CO 


"U 

c 

CO 

0 ) 


W ■§ 

o Q 

*5 

<s O 

3 4 ! 

o at 
o c □ 

CO C t' 

a o ^ 

(0 o £ 


* 

C 

o 

sa >> 

> J 5 
_c C/3 
<0 > 
o > 

* “5 CD 

V-o £ 
t c -K 
0) as ro 
.c 

TO © T 3 

3 O) 2 

° 0 m 

08 ? 

£ Q ® 
co c 0 
co jo □ 

3< 

CO ^ 
CD 0) 
^ CD 

!t= Q " 
0 


CD 

O 

o 

CM 

</> 


Q. 

O 

CO 


-T' v 

« § 


CD 

co 2 

■— o 
<0 CM 

CO 

s < £ 


CO 


CO 

3 

D) 


CO 

"co 

E 

u. 

a> 


CO 

3 

c 

c 

< 

JC 

r^ 


_G) 

'> 

J 2 

c 

3 


o 

LL 

CO 

2 

< 

CO 

< 

z 


co 

, ^r 


03 

uj 

o 

LL 

CO 


CD 

O 

O 

CN 

CO 


A Technical Fellow, ERCJnc., MSFC, Huntsville, AL 
! CFDRC, Inc., Huntsville, AL 



£ C/> 

o '</> 

!E >* 

0 > 05 
> C 

<D < 

3 t 

3 O 
s: cl 
CO £ 

to a 


CD 

> 

■ ■■■ 

O 

CD 

■ ■■ 

-Q 

o 


o .<2 

0 ^ 
"D 0 

o o 
E st 



co £ 

I « 

«“■ 

o| 

0 3 
> 0 
0 —I 

O £ 

I— 0 




CD 

1.1 
*— -I—* 

^ o 
c 

c 
o 
o 


o c 

0 
sz 


I- c 

(!) 

+■» 

O o 


0 

o 

t o 
0 +- 
Q_ 0 

■gE 

■Q O 

0 t 
■O o 

0 Q- 
0 

ii 


CO 

'i— 

-Q 

0 

"O 

0 


0 

O 

3 

“O 

O 


c 

0 

E 

c 

o 


> 
c 

Q. 0 



(N 


CD 

O 

o 

CNI 

0 


Vehicle Debris Environment 

(Impact kinetic energies, velocities, angles, locations) 



A Full 3-D CFD Model 



it 

0 

1 

£ 
■ MB 





+ 





OL 

a: 



03 

(/) 


LU 

§ 

(0 

CO 


CO 

O' 

LU 


O 

LL 

c o 


CD 

O 

O 

CM 

CO 



0 

'i_ 

CD 

0 

0 

i_ 

0 

> 

o 


o 

CZ 

0 

_l 

E 

o 

s_ 

M— 

L_ 

0 

> 

o 

0 

Q 

LL 

O 


o 

t 

0 

> 

O 

0 

i- 


</> 
«*— » 

c 

o 

Q. 

’l. 

O) 

c 

o 


o 

I s - 


Q. 

0 

JO 

O) 

“D 

0 

v_ 

3 

O 

3 

-*— • 
(/) 

"D 

0 


"D 

O 

-Q 


O 

CO 

0 

-C 

-I— ' 

“D 

c 

05 

O 

LL 

CO 


05 

0 

0 

0 

_3 

O 

O 

CL 

«4— 

O 

c 

o 

to O 
.£ CL 

F < 

i« 

Jf to" 

g l 

|o 

o >< 
0 ~ 
X 5 
0 < 

0 

c 

o 


o 

0 


(/) 

c 

o 

05 

0 

o 

o 

o 

cT 

CM 

A 


0 

O 

c 

0 

CD 

i— 

0 

> 

c 

o 

o 


c 

3 


CD 

C 

'c 

o 

U— » 

d 

c 

o 

0 

1 

0 

I— 

CL 

O) 

c 

CO 


o 

d 

ii 


c 

£ 

o 

"O 

i— 

0 

_Q 

E 

3 

C 

-C 

o 

05 


05 


C 

3 


0 

-Q 

0 

“O 

05 

■D 

C 

CO 

0 

■D 

O 

O 

tr 

o 

CL 0 
CO “O 

c o 

05 O 

-t 0 
0 o 

O CO 

tg fe 
00 -S 
0-0 

0 

■3 0 
_Q CD 
0 05 

■o E 

05 05 
r- "O 


c 

o 


o 

CO 

Ql 

E 


0 

-Q 

05 


O 0 

§ o 

E* ^ 

c 

to 


o 

LL 

— CO 
0 ^ 
■D ^ 
O > 

E 

O ? 

< CD 

O 0 

= >§ 
o 13 
0 0 

^ Q 

O ^7 

0 o 
0 

E o 

2 "D 

— 0 

0 tr 


0 

> 

o 

o 


o 


"O 

o 
E 

Q 
lj- O 
O 0 
0 ^ 

.2 E 
j= o 
0 4 = 
> 


0 
0 
■o 
o 
o 
-*— » 
o 
0 
Q_ 

I| 

‘I® 

*= o 


co-g 

0 ^ 
O 0 

E § 

0 s_ 

-go 

2 co 

0 D 

ro -o 

0 C 

o 0 

'€ O 

s°e 

■= £ 
-Q O 
0 4 = 

Q 


O 

LL 

0 

>% 

_Q 

O 

'-I— * 

0 

i_ 

CD 

0 

C 

St 

O 


O 

LL 

0 

E 

o 


0 

0 




O 

CN 

0 


* Overflow code by Pieter Buning, LaRC 


Shuttle 3-D Exhaust Plumes 



</) 

CD 


E 


3 

Q. 



o 

-o +- 

O 0) 
o « 

% O -a 
o > c 
C O D1 

o 




CD 

* 

CO 


■O 

O 

o 


■O 

CO 

■O 


co 

+-» 

C/> 

LL 

< 



0 25 50 75 


(!33d)30NVlSia 1VIXV 


a> 


r 

2 


o 

CD 



(/) 



H 





Pre-Launch before SSME Start 







SSMEs Started and at 100 % Power Level 
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Uncertainties being analyzed 

Cycle 2 implementation and uncertainty quantification to begin 
September 2006. 
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